Optomechanical coupling in a two-dimensional photonic crystal defect cavity 
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Periodically structured materials can sustain both optical and mechanical modes. Here we inves- 
tigate and observe experimentally the optomechanical properties of a conventional two-dimensional 
suspended photonic crystal defect cavity with a mode volume of ~3 (A/n) 3 . Two families of mechani- 
cal modes are observed: flexural modes, associated to the motion of the whole suspended membrane, 
and localized modes with frequencies in the GHz regime corresponding to localized phonons in the 
optical defect cavity of diffraction-limited size. We demonstrate direct measurements of the optome- 
chanical vacuum coupling rate using a frequency calibration technique. The highest measured values 
exceed 250 kHz, demonstrating strong coupling of optical and mechanical modes in such structures. 

PACS numbers: 42.70.Qs, 43.40.Dx 



Cavity optomechanics [1] exploits the coupling of me- 
chanical oscillators to the light field via radiation pres- 
sure. On the applied side, such coupling may be used to 
enable novel radiation pressure driven clocks [2], make 
highly sensitive displacement sensors, tunable optical fil- 
ters [3], delay lines [4] or enable photon storage on a 
chip [5, 6]. On a fundamental level, such systems can 
be exploited for demonstrating that nano- and microme- 
chanical oscillators can exhibit quantum mechanical be- 
havior [7]. It has been shown experimentally that it 
is possible to cool a mechanical oscillator intrinsically 
via radiation pressure dynamical backaction [8]. In or- 
der to reach the ground state of mechanical motion and 
enable manipulation in the quantum regime, one ap- 
proach consists in down-sizing the oscillator, thus shift- 
ing the quantum-classical transition towards higher tem- 
peratures. Among various sub-micron optomechanical 
systems presently investigated [9] , suspended membranes 
containing a photonic crystal cavity offer strong light con- 
finement in diffraction-limited volumes and are therefore 
natural candidates for achieving strong optomechanical 
coupling. Recently, optomechanical coupling in ID pho- 
tonic crystal systems [10] has been observed in patterned 
single [11] and dual nanobeams (zipper cavities) [12, 13]. 
It would be highly desirable to extend such optomechan- 
ical coupling to 2D systems, notably photonic crystal 
defect cavities. Such cavities offer the strong light con- 
finement possible (i.e. small mode volume), high quality 
factor (Q) [14] and have been used for studying cavity 
QED using quantum dots [15] or for realizing nanolasers 
[16, 17]. Recently, a 2D optomechanical photonic crys- 
tal slot cavity has been reported [18]. While mechanical 
displacement due to strong radiation force generated by 
band-edge modes in bilayer photonic crystal slabs has 
been reported as well [19], to date optomechanical cou- 
pling of the 2D conventional defect cavity has not been 
studied. 



In this paper, we demonstrate optomechanical coupling 
using a photonic crystal defect (L3) cavity. We provide a 
direct and robust experimental determination of the vac- 
uum optomechanical coupling rate [30] using frequency 
modulation, showing a particularly strong coupling for 
the localized mechanical modes, which may also be cou- 
pled to quantum dots in future studies [20] . 




FIG. 1. (a) Scanning Electron Microscope side view of the 
cavity, (b) Microphotoluminescence spectrum of the pho- 
tonic crystal slab cavity obtained under non-resonant contin- 
uous optical excitation at normal incidence with an excitation 
power of 100 fiW at 532 nm (c) Micrograph (false colors) of a 
defect cavity fiber-taper system used to read out mechanical 
motion of the cavity, (d) Experimental setup (ECDL: Exter- 
nal cavity diode laser, OI: Optical Isolator, EOM: Electro- 
optical modulator, FPC: Fiber polarization controller, PhC: 
Photonic crystal defect cavity, PD: Photo diode). 

The optomechanical device under study consists of 
a 262-nm thick InP suspended membrane containing a 
two-dimensional photonic crystal defect cavity shown in 
Figure 1(a). The cavity, following the design proposed 
in [14], contains three missing holes in a line of a perfect 
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triangular lattice of holes with a lattice constant of 
a = 430 nm and a radius r — 90 nm. At both edges of 
the cavity the holes are displaced outwards by d = 0.18a, 
in order to obtain high optical quality factors. The 
cavity is fabricated using electron beam lithography, 
inductively coupled plasma etching [21], and wet etching. 
The cavity incorporates a single layer of self-assembled 
InAsP quantum dots [22] at its vertical center plane for 
cavity characterization. The whole structure is grown by 
metalo-organic chemical vapor deposition. The quantum 
dot density is ~15xl0 9 cm -2 , and their spontaneous 
emission is centered around 1560 nm at 300 K with 
an inhomogeneous broadening of about 150 nm. The 
presence of the dots inside the cavity allows to identify 
the spectral properties of the fundamental optical mode 
of the cavity by photoluminescence measurements [23], 
as shown in Figure 1(b). The resonance wavelength of 
the fundamental mode is centered around 1555 nm and 
the cold-cavity quality factor is measured to be ^10 4 
(cavity linewidth is k/2it w 20 GHz). The suspended 
photonic crystal membrane lies on top of a 10 /im high 
mesa structure (see Fig. 1(a)). The mesa structure is 
processed to enable positioning of a tapered optical fiber 
in the evanescent field of the cavity, while precluding any 
interaction with the nearby substrate. The shape of the 
membrane resembles a Bezier curve, which was chosen to 
increase optomechanical coupling to the flexural modes. 

The setup used in the experiment is depicted in Figure 
1 (d) . An external-cavity diode laser is used for the read- 
out of the mechanical motion. Coupling to the optical 
modes of the suspended membrane is achieved with the 
optical fiber-taper technique [24] . Piezoelectric actuators 
enable an accurate positioning allowing to optimize the 
gap between the fiber-taper and the defect cavity, and 
thus to increase evanescent coupling. Despite careful re- 
duction of the taper-cavity gap, only a small fraction of 
the light can be coupled into the defect cavity, typically 
not exceeding 10 % of the incoming laser power. The 
cavity-fiber system is kept in a vacuum chamber with a 
pressure below 1 mbar. Laser light coupled inside the 
photonic crystal cavity leads to local heating, which in- 
duces a thermal effect arising from the temperature de- 
pendent refractive index n [25]. As dn/dT > (T is the 
temperature) for InP [26] , the region detuned to the blue 
side of the resonance allows thermal passive locking [25] . 
In our experiments the laser frequency is chosen to corre- 
spond to the blue-detuned side of the fringe of the optical 
mode requiring no further locking. Mechanical motion of 
the membrane is imprinted on the transmitted optical 
intensity through modulation of the internal cavity field. 
An electro-optical modulator is used for frequency mod- 
ulation to determine the optomechanical coupling rate. 
The transmitted signal is detected by a fast receiver and 
the electrical signal is analyzed with an oscilloscope as 
well as an electronic spectrum analyzer, which is used 



for the spectral analysis. 
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FIG. 2. (a) Detected frequency noise spectrum in the 1 MHz 
- 200 MHz range presenting a series of peaks corresponding 
to the different mechanical modes labeled by numbers (Black 
curve). The red curve represents a spectrum acquired with 
the laser being detuned out of resonance. Inset: Calibrated 
frequency noise spectrum of the fundamental mode (#1) with 
a Lorentzian fit (red line), (b) Spatial displacement pattern 
of the first eight mechanical modes and the prominent mode 
around 150 MHz, as obtained from finite element modeling. 

For a launched laser power of 1.3 mW more than 20 
mechanical modes are observed in the frequency range 
between 10 MHz and 1 GHz. These modes can be sepa- 
rated into two mode families. The first family consists of 
flexural modes present in the low-frequency range (below 
200 MHz), whereas the second family consists of local- 
ized modes. Flexural modes, whose spectrum is shown 
in Figure 2(a), correspond to the movement of the whole 
membrane. In order to identify the various modes, we 
modeled the mechanical properties of the photonic crys- 
tal slab structure by finite element modeling (COMSOL 
Multiphysics). Realistic geometry parameters were taken 
into account, including the under-etching of the mesa 
structures between which the membrane is suspended. 
A good agreement between measurements and modeling 
is obtained using a Young modulus of 20 GPa (slightly 
smaller than usual values observed in bulk InP materials 
[27] attributed to the perforation for the photonic crys- 
tal). Figure 2(b) shows the displacement patterns of the 
first eight modes as well as a prominent mode around 150 
MHz. 

Localized modes, shown in Figure 3, correspond to me- 
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FIG. 3. (a) Detected frequency noise spectrum in the 200 
MHz - 1.1 GHz range presenting a series of peaks correspond- 
ing to the different mechanical modes labeled by numbers 
(Black curve). The red curve represents a spectrum acquired 
with the laser being detuned out of resonance. CP denotes 
the calibration peak resulting from the frequency modulation 
measurement. Inset: Measured versus simulated frequency 
for mechanical modes which could be assigned (#1-15 and 
#19-21). (b) Spatial displacement pattern for the first four 
orders of localized mechanical modes, as obtained from finite 
element modeling, (c) Simulated (blue) and experimentally 
determined (red) progression of the resonance frequency of 
the localized modes versus mode number (mode #22 was de- 
termined in a separate measurement). 



chanical displacement of the membrane localized in the 
cavity core of the defect. We were able to resolve the 
fundamental localized mode at 0.46 GHz as well as the 
three higher mode orders (at 0.72 GHz, 0.99 GHz and 
1.26 GHz). The inset of Figure 3(a) compares the simu- 
lated against the measured resonance frequencies of the 
mechanical modes that could be assigned, revealing ex- 
cellent agreement both for flexural and localized modes. 
The progression of the resonance frequency of the local- 
ized modes versus mode number is shown in Figure 3(c) 
and follows a linear behavior with mode number. 

Importantly, localized mechanical modes coincide spa- 
tially fully with the optical defect cavity mode. There- 
fore, the photonic crystal not only offers strong optical 
confinement, but simultaneously ultra-high phonon con- 
finement. We note that the localization occurs in the 
absence of a phononic band gap [28]. Due to the co- 
localization of the optical and mechanical mode within 
the defect cavity, strong optomechanical coupling is ex- 
pected. This is visualized in Figure 4(a), where the re- 
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FIG. 4. (a) Simulated distribution of the electromagnetic 
field (left), Scanning Electron Microscope image of the de- 
fect (middle) and simulated third order localized mechanical 
mode (right), (b) Calibrated frequency noise spectrum of the 
third order localized mechanical mode #21 (blue points) with 
a Lorentzian fit (red line). A vacuum optomechanical rate of 
go/2n = 268 kHz is determined. 



suits of the mechanical FEM simulation for the third or- 
der localized mode are compared with the spatial dis- 
tribution of the electromagnetic field of the optical de- 
fect mode, which was obtained through a finite differ- 
ence time-domain (FDTD) simulation. The high overlap 
of the mechanical displacement with the distribution of 
electromagnetic energy promises strong optomechanical 
coupling. In a second experiment wc determined the cou- 
pling rate for various flexural modes and the second and 
third order localized modes. 

Usually, the optomechanical coupling strength is de- 
termined by two parameters - the optomechanical cou- 
pling parameter G — ^f, with uj c being the resonance 
frequency of the optical resonator and x denoting the dis- 
placement of the mechanical oscillator, and the effective 
mass m c g [29] of the mechanical mode. The necessity of 
introducing the effective mass routinely arises from the 
arbitrary definition of x, which often cannot be consis- 
tent with the displacement pattern of different mechani- 
cal modes of the system. Particularly for photonic crys- 
tals it is difficult to define an unambiguous displacement 
direction of a mechanical mode. One attempt to circum- 
vent this problem is to introduce an effective length L e g 
determined by a perturbative expression of the overlap 
of mechanical displacement and the electromagnetic field 
distribution [12, 13]. Recently, it was suggested that the 
vacuum optomechanical coupling rate go would be a more 
proper quantity for optomechanical systems [30, 31]. In 
analogy to cavity Quantum Electrodynamics (cQED), go 
is defined as go — G ■ x zp i, with x zp f = \J h/2m c gfl m be- 
ing the zero-point-fluctuations of the mechanical oscilla- 
tor (h indicates the reduced Planck constant and Q m /2ir 
the mechanical resonance frequency) . As all relevant op- 
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tomechanical parameters can be derived through knowl- 
edge of go > acquiring its value would make the determi- 
nation of G and m e s redundant. As described in [30], 
the value g can be determined experimentally via 

.90 = \jSujoj (^m) ^T 1 , (1) 

with n being the average phonon occupancy of the 
mechanical mode, T m /2ir being the mechanical damping 
rate and S uu (fl m ) being the frequency spectral density 
of cavity frequency noise evaluated at the mechanical res- 
onance frequency. For high phonon occupancy one can 
approximate n ~ k-QT/Ml m >• 1, with as Boltzmann's 
constant. S uu (fi TO ) can be experimentally determined by 
a frequency modulation technique [32] . 

We performed a measurement of the cavity frequency 
noise by applying a known phase modulation to the laser 
using an LiNbC>3 electro-optical modulator (cf. Figure 
1). The incoming laser power was reduced to 0.9 mW. 
To make sure that the mechanical modes were not driven 
thermally, the laser frequency was slightly detuned from 
the side of the fringe in both directions just before the 
experiment. With no change in the resonance frequency 
of the mechanical mode occurring, we inferred that the 
mechanical modes were only driven by thermal Brown- 
ian motion (T=300 K). The modulation frequency was 
chosen to be close to the resonance frequency of the me- 
chanical mode to be calibrated. The measurement also 
allows calibration of the frequency noise produced by the 
mechanical mode. A detailed account of the calibration 
method is given in [33]. A calibrated frequency noise 
spectrum for the fundamental flexural mode is shown in 
the inset of Figure 2(a), and the optomechanical vacuum 
coupling rate was determined to be go/2ir — 234 Hz. 
The mechanical quality factor for this mode is Q m = 890 
being the highest for all flexural modes. The coupling 
rate for the flexural modes increases with the mode num- 
ber up to several kHz. The optomechanical coupling for 
the localized mechanical modes was determined to be 
g /2ir = 199 kHz (Q m = 160) and g /2ir = 268 kHz 
(Qm = 180) for the second and third order, respectively. 
A calibrated frequency noise spectrum for third order lo- 
calized mode is shown in Figure 4(b). The high values of 
go/2n give a definite experimental proof of the high op- 
tomechanical coupling between a photonic crystal defect 
cavity and a localized mechanical mode. These coupling 
values are two orders of magnitudes higher than mea- 
sured in whispering gallery mode toroidal resonators [34] 
and doubly-clamped strained silicon nitride beams in the 
near-field of a silica toroidal resonator [9] , both of which 
are ca. go/2ir s=s 1 kHz. Moreover, the measured values 
are as high as the recently reported coupling of a flexural 
mode to a photonic crystal slot cavity [4]. The values of 
go/2-K which could be unambiguously obtained are spec- 
ified in Table I. 



TABLE I. Optomechanical vacuum coupling rates for various 
modes determined experimentally via the frequency modula- 
tion technique. 



Mode 
index 


Measured frequency 
n m /2?r (MHz) 


Measured vacuum 
coupling rate go/2-K (kHz) 


1 


9.54 


0.23 


2 


24.45 


0.67 


5 


46.87 


2.26 


8 


59.83 


7.26 


20 


716 


199 


21 


991 


268 



The calibration technique used has the clear advantage 
that it allows to determine go for any mechanical mode 
from the mere knowledge of its mechanical linewidth and 
the mechanical mode occupancy. The method does not 
require knowledge about the optical linewidth, the spe- 
cific transduction mechanism of the signal through the 
optical cavity, the spatial distribution of electromagnetic 
energy or the displacement pattern of the mechanical 
mode. Thus the method is particularly suitable for 
optomechanical calibration in photonic crystals due 
to the frequently complex spatial distribution of both 
mechanical and optical modes. Moreover, the present 
planar architecture can be used for photon-phonon 
conversion experiments [4]. 

In conclusion, we demonstrated optomechanical cou- 
pling in a two-dimensional III-V photonic crystal defect 
cavity. We observed both flexural as well as localized 
mechanical modes. Furthermore, we provide direct mea- 
surements of the vacuum optomechanical coupling rate 
in a photonic crystal and measure for the first time the 
coupling of a localized photonic crystal defect mode with 
a localized 2D mechanical mode. Coupling rates for the 
localized modes exceed 250 kHz and arc two orders of 
magnitude larger than in conventional optomechanical 
systems. By integrating a single quantum dot in the de- 
fect cavity, a variety of experiments can be envisioned 
such as laser cooling [20] and coupling of a quantum me- 
chanical oscillator to an artificial atom, once the ground 
state of the mechanical oscillator is reached. 

Funding for this work was provided by European 
NanoSci-ERA project NanoEPR, through QNEMS and 
MINOS by the FP7, the NCCR Quantum Photonics, 
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Supplementary information - Optomechanical coupling in a two-dimensional photonic 

crystal defect cavity 



DETERMINATION OF THE OPTOMECHANICAL VACUUM COUPLING RATE 



The optomechanical vacuum coupling rate go/2ir is determined via a frequency modulation technique. An electro- 
optical modulator (EOM) is used to phase modulate the laser carrier before its coupling into the cavity. The modulated 
phase $ of the signal can be written as 



=co t + l3cos{Cl mod t), (SI) 

with u as the radial laser frequency t as the time, (3 as the phase shift factor and f2 mo( z/27r as the modulation 
frequency. The transduced signal from the photonic crystal cavity is detected with a fast receiver and the spectrum is 
resolved with an electrical spectrum analyzer (ES A) . The spectrum exhibits the Lorentzian spectra of the mechanical 
modes and a calibration peak at tt mod /2Tr resulting from the modulation [1, 2]. If the modulation frequency is chosen 
to be close (a couple of mechanical linewidths) to the resonance frequency f2 m /27r of the mechanical mode to be 
calibrated, go can be determined as [3] 

2 ^ \ P 2 £lmod r m (Ct m ) 

% ~ 2n 2 4 • RBW S uul (n mod ) ' 

with n as the average phonon occupancy, r m /27r as the dissipation rate of the mechanical oscillator, RBW as the 
resolution bandwidth of the ESA, 5 WW (f2 m ) as the double-sided frequency spectral noise density evaluated at the 
mechanical resonance frequency and S uu (ti mod ) as the double-sided frequency noise spectral density evaluated at the 
modulation frequency. Only the proportion between S uu (Sl m ) and (Jl m od) is required to be known, and it can 
be readily obtained from the measured peak power spectral density values of the mechanical and calibration signals. 
The power spectral density is usually given in dBm in the spectrum of the ESA. The proportion can be found with 
the following relationship 

s uu (On) = s uu (n mod ) • i -( p — *- p -)/ 10 , (s 3) 

with P rno d being the power at the peak of the calibration signal expressed in dBm and P m being the power at the 
peak of the mechanical mode, which is expressed in dBm as well. 



(S2) 



The frequency modulation technique can be also used to calibrate the spectrum in absolute units. This is done by 
directly determining S uu (fl mo d) via [3] 

(timed) = ( f .^b W • ( S4 ) 

By using an equivalent scheme depicted in Equation S 3, the spectrum can be calibrated. The calibration holds true 
as long as the frequencies of the values to be calibrated are close enough to the resonance frequency of the calibration 
peak. Otherwise, a generalized transduction coefficient needs to be included [3]. 



CALCULATION OF CORRECTION FACTORS FOR THE ESA SPECTRUM 



When measuring electrical signals with an ESA, one usually expresses the power detected at a certain frequency in 
the logarithmic units of dBm. By averaging one obtains an average of the logarithmic power, which is not equal to the 
logarithmic expression of the averaged power - the value that is of significance in the measurements. The discrepancy 
is different for random signal, such as background noise and the mechanical spectrum (resulting from thermal motion), 
and for coherent signals, such as a calibration peak with a large enough signal-to-noise ratio. For random signals it 
is known [4] that the measured value is about 2.5 dB below the actual one, whereas for truly coherent sources the 
measured signal corresponds to the real one. For signals that have contributions both from a coherent source and 
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random signals, such as a calibration peak with a rather low signal-to-noise ratio, the correction factor depends on 
the magnitude of both contributions. 

Modern ESAs can often compensate this discrepancy internally, however as this function was not applied during 
our measurements, we need to use the correction described above. In the following we will derive the general 
expression for the correction of the measured signal. 

A signal arriving on the ESA can be decomposed in its in-phase (/) and out-of-phase (Q) part. Spectrum analyzers 
respond to the magnitude of the signal within their RBW passband [4]. The magnitude of a signal (voltage) v 
represented by / and Q is given by 

v= Vl 2 + Q 2 - (S5) 
The average power arriving on the detector is given by 

where an input impedance of 50 f2 is considered (equations are generally dimensionless in this section) and the 
brackets () imply that the mean value is taken of the expression inside them. In the average mode the ESA averages 
the logarithmic input signal, thus leading to the following expression for the measured power 

Pmeas = (lO-\og(^r^^y (S 7) 

The logarithmic value of the actual power of the signal is given by 

P log = 10 • log (P) (S8) 

It should be intuitive that the values for P log and P meas are not equal, just as the log of the average is not the same 
as the average of the log. 



The mean values given in Equations S 6 and S 7 can be determined by multiplying the variable being measured 
with its probability density function (PDF) and integrating over the possible values of the variable. The PDFs of the 
quadratures / and Q can be assumed to be Gaussian with a certain variance 07 for I and oq for Q, and a certain 
offset of the distribution from being /ij for I and /j,q for Q. In general, one can assume the variance to be the same 
for both quadratures, so that oj — (Jq — a. For random signals it holds that f-ti = Hq = 0, whereas for a coherent 
signal one has [1q = 0, but fij = [i > 0. 



Using the assumptions stated above, one can write 



PDF(,) = d^ exp (-W (S9) 



and 



PDF (Q) = —7=- cxp I ) • (S10J 



From the Equations S6 and S 7 it is obvious that we are interested in the mean values of v 2 . This requires one 
to find the PDFs for this variable from the PDFs of / and Q. This is accomplished by using a general relationship, 
which links the PDFs fx t (xi) of independent random variables Xj, i = 1,2, ...n to the PDF /y (y) of some variable 
Y = G (X\, X2, ■■■X n ). This relationship is given by 

/OO /'OO p 00 

/ ... / dx 1 dx 2 ...dx n fx 1 { x i)fx 2 {x2)-fx n (xn)S(t-G(x 1 ,x 2 ,...x n )), (S 11) 

-00 J —00 J -co 



where S denotes Dirac's delta function. 
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Using Equation S 11 one can determine the PDF of v 2 to be 

™m-£«p(-^)-*GS). 

with J as a Bcssel function of the first kind. 

Using this expressions for the PDF of v 2 , one can calculate P and P meas by using Equations S 6 and S 7. For P 
one obtains 

v 2 

■ ■ PDF (v 2 ) 

50 ' < S13 > 
and for Pmeas one obtains 



p 2 + 2a 2 



-Pmeas (M, ^) = ( 10 • log 



50 



Jo 



dv 10 • log ( ^ j • PDF (w 2 ) 



io -(r (o,^)+in(g 

In (10) 

where T ^0, ^ ) is an incomplete Gamma function and In denotes the natural logarithr 
For random signals we set fi to obtaining 



iiu 10 -( r (°'^)+ ln (^)) 

io(r B + in(||)) 



Pmeas {P- = 0, Cr) = Hm P mea s (M, O") = lim 

ai->o ^-s-o m(10) 



(S14) 



(S15) 



In (10) 

with as the Euler-Mascheroni constant (r^ w 0.577216). 

Finally, we obtain a general expression for the discrepancy between the actual and the measured signal 

A ( M , a) = P roeos _ 10 . log (P) = V V 2CT J (1Q) ^+ 2 ^ . (s 16) 

For the special case of a random signal the final expression reads 

A( M = ) . ) = limA( M ,. ) = -^, (S17) 

The discrepancy in this case is a constant value being approximately A (p, = 0, a) = —2.50682 and independent of 
a as expected. 

For a strongly coherent signal with \i ^> a one obtains A (/i 3> cr, cr) = 0. This means that for a strong coherent 
drive there is no discrepancy between the measured and the real one. 
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FIG. SI. Reflection coefficient Teom of the power being received by the combined EOM+cable system versus frequency 
measured with a network analyzer. 
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FIG. S2. Transmission coefficient Tbnc of the power being transmitted through the BNC cable linking the signal generator 
and the EOM versus frequency. The measurement was performed with a network analyzer. 

EXPERIMENTAL IMPLEMENTATION 



The phase shift factor j3 can be determined by 

P = V/V« ■ it = ^WfZir/V^ (S 18) 

with V as the applied voltage, V n as the voltage needed to induce a phase shift of 7r, P r j as the input power of 
the signal generator providing the modulation signal and Z as the impedance between the signal generator and the 
EOM. To have a quantitavie understanding of the actual voltage applied on the EOM, we directly measured the 
reflectance of the system consisting of the EOM and the BNC cable linking it to the signal generator. Moreover, we 
measured transmittance through the BNC cable to account for potential losses. The measurements were performed 
using a network analyzer. Figure SI shows the fraction of the reflected power Teom of the combined EOM+cable 
system as a function of the output frequency of the signal generator. Figure S2 shows the fraction of transmitted 
power T BNC through cable as a function of frequency. It is obvious that for high frequencies the power loss becomes 
fairly strong, which is due to radiation losses in the BNC cable. This fact is undermined by Figure S2 showing the 
reflection coefficient Tbnc through the cable, which is well below the one of the total EOM+cable system. 
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FIG. S3. Reflection coefficient Fbnc of the power being transmitted through the BNC cable linking the signal generator and 
the EOM versus frequency. The measurement was performed with a network analyzer. 

The factor f3 can be thus determined as 

P = ^2P rf • (1 - 10 r *Wio) . ioiw/io . Z out Tr/V^. (S 19) 

Measurements gave a value of V„ « 7 V which is consistent with the value provided of the manufacturer of the 
EOM. 

A spectrum of the fundamental flcxural mode is shown in Figure S4 together with the calibration peak. As discussed 
in the previous section, the measured spectrum needs to be readjusted. This is done by first fitting the mechanical 
with a Lorentzian and the calibration peak with a Gaussian. The noise level is obtained from the baseline of the 
Lorentzian, and this value is subsequently used to determine the variance a as defined in Equations S 9 and S 10 
through use of Equation S 15. From the Gaussian fit we obtain the maximum value of the calibration peak. Using 
this value together with the determined value of a, one can numerically solve Equation S 14 to obtain /i. The value 
of /i 2 /50 is the corrected maximum value of the calibration peak. The spectrum of the random signals, consisting 
of the background noise and the mechanical response, is corrected with the constant value given in Equation S 17. 
The corrected values for the maxima of the mechanical spectrum and the calibration peak are used together with 
r m , inferred from the Lorentzian fit, to determine go as described in Equation S2. By using Equation S4 one can 
calibrate the spectrum in absolute frequency units. 



The values of go/^n, that could be unambiguously determined, are summarized in Table SI. The value of go/2ir 
without the correction for logarithmic averaging is also stated together with the calculated correction for the calibration 
peak. The required correction for the calibration peak is fairly small for most measurements due to the large enough 
signal-to-noise ratio. In those cases the main correction stemmed from the adjustment of the random signals by a 
constant factor of around 2.5 dBm. 



Mode 


Q.m/2-K 


g /2n 


go/2-K (without correction) 


Correction factor 
for the calibration 
peak 


1 


9.54 MHz 


234 Hz 


175 Hz 


0.000 dBm 


2 


24.48 MHz 


668 Hz 


500 Hz 


0.004 dBm 


3 


46.87 MHz 


2.26 kHz 


1.69 kHz 


0.007 dBm 


4 


59.83 MHz 


7.26 kHz 


5.44 kHz 


0.000 dBm 


5 


716.39 MHz 


199.1 kHz 


148.4 kHz 


0.045 dBm 


6 


990.88 MHz 


267.8 kHz 


193.4 kHz 


0.322 dBm 



TABLE SI. Table with go for different mechanical modes as well as different calculation conditions. The last column specifies 
the correction factor for the peak height of the calibration peak. 
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FIG. S4. Mechanical mode (blue) with a calibration peak (green) including the respective fits. 



The values of go/27r, that could be unambiguously determined, are summarized in Table SI. The value of go/2ir 
without the correction for logarithmic averaging is also stated together with the calculated correction for the calibration 
peak. The required correction for the calibration peak is fairly small for most measurements due to the large enough 
signal-to-noise ratio. In those cases the main correction stemmed from the adjustment of the random signals by a 
constant factor of around 2.5 dBm. 
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